Large-scale growth evolution in the Szekeres inhomogeneous cosmological models with 

comparison to growth data 



Austin Peel*, Mustapha Ishak*, and M. A. Troxel* 
Department of Physics, The University of Texas at Dallas, Richardson, Texas 75083, USA 

(Dated: December 12, 2012) 

We use the Szekeres inhomogeneous cosmological models to study the growth of large-scale struc- 
ture in the universe including nonzero spatial curvature and a cosmological constant. In particular, 
we use the Goode and Wainwright formulation of the solution, as in this form the models can be con- 
sidered to represent exact nonlinear perturbations of an averaged background. We identify a density 
contrast in both classes I and II of the models, for which we derive growth evolution equations. By 
including A, the time evolution of the density contrast as well as kinematic quantities of interest can 
be tracked through the matter- and A-dominated cosmic eras up to the present and into the future. 
In class I, we consider a localized cosmic structure representing an overdensity neighboring a central 
void, surrounded by an almost Friedmann-Lemaitre-Robertson- Walker background, while for class 
II, the exact perturbations exist globally. In various models of class I and class II, the growth rate 
is found to be stronger in the matter-dominated era than that of the standard lambda-cold dark 
matter (ACDM) cosmology, and it is suppressed at later times due to the presence of the cosmo- 
logical constant. We find that there are Szekeres models able to provide a growth history similar 
to that of ACDM while requiring less matter content and nonzero spatial curvature, which speaks 
to the importance of including the effects of large-scale inhomogeneities in analyzing the growth 
of large-scale structure. Using data for the growth factor / from redshift space distortions and 
the Lyman-a forest, we obtain best fit parameters for class II models and compare their ability to 
match observations with ACDM. We find that there is negligible difference between best fit Szekeres 
models with no priors and those for ACDM, both including and excluding Lyman-a data. We also 
find that the standard growth index 7 parametrization cannot be applied in a simple way to the 
growth in Szekeres models, so a direct comparison of the function / to the data is performed. We 
conclude that the Szekeres models can provide an exact framework for the analysis of large-scale 
growth data that includes inhomogeneities and allows for different interpretations of observations. 

PACS numbers: 98.80.Es,98.80.-k,95.30.Sf 



I. INTRODUCTION 



Studying inhomogeneous cosmological models is becoming increasingly important as the available cosmological data 
expand and improve. These exact solutions to Einstein's field equations provide frameworks to analyze the data with 
a wider range of possible interpretations. They can also be used to represent nonlinearities, which cannot be described 
using a Friedmann-Lemaitre-Robertson- Walker (FLRW) model plus linear perturbations. The real universe is lumpy 
and exhibits a striking variety of small- and large-scale inhomogeneous structures, such as large voids, clusters, and 
superclusters of galaxies. Some of these superstructures can be as large as 5% or 10% of the Hubble scale, including 
the well-known Pisces-Cetus Supercluster Complex and the Sloan Digital Sky Survey Great Wall. 

Exact solutions to Einstein's equations representing inhomogeneous cosmological models have been the subject of 
a number of theoretical studies [lj, |2|; however, their comparison to observational data is only a recently emerging 
field. For example, although spherically symmetric, the Lemaitre-Tolman-Bondi (LTB) solutions have achieved some 
proof-of-principle successes in their application to cosmology, where certain models have been shown to be able to fit 
type la supernova and baryon acoustic oscillations (BAO) data, as well as the cosmic microwave background power 
spectrum pMLa]. 

We use in this work the Szekeres inhomogeneous cosmological models. These exact solutions have no Killing vectors 
(i.e., no symmetries) and are well suited to represent the lumpy universe we observe. They have a broad potential in 
cosmology to provide a more realistic description of the Universe 0, [TtJ ■ The solution was first derived by Szekeres 
pH . [l9j from a general metric with irrotational dust as the source of the spacetime. Several authors have studied the 
models analytically and numerically; a partial list includes [ljl [2fj| - l36j . The Szekeres models have been extended to 
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include pressure in (37j but with some known limitations on the state variables, and they were also extended by (38j 
to include a cosmological constant. 

The expansion and growth histories must both be consistent with observational data in any cosmological model 
that purports to describe the universe. Cosmological distances are necessar y to understand the expansion history, and 
distances versus redshift have been previously studied for Szekeres models [32l l34l 39] . The growth rate of large-scale 
structure in the matter-dominated era for flat and curved cases of class I and class II Szekeres models have been studied 
in our earlier paper [2(|. In that work, a formulation of the models due to Goode and Wainwright (GW) was used 
[22, [23], which is well suited to studying structure growth in the universe. In the GW representation, Szekeres models 
can be interpreted as exact nonlinear perturbations of some smooth associated FLRW background. Recently, the 
authors of [35| extended the work of GW for class II flat models with a cosmological constant and examined nonlinear 
inhomogeneities in a ACDM background. Specializing to the spatially flat case, they showed that the inclusion of a 
cosmological constant allows models in some cases to avoid shell-crossing singularities. The evolution of structure in 
Szekeres models has also been studied via invariant density contrast indicators [2l| [2l| [3(| , but we take a different 
approach here. 

The present paper generalizes our previous work [20I ] to include a cosmological constant in exploring the growth 
history of Szekeres models. Thus, the analysis now applies fully to the time since recombination up to the present and 
into the future, both dominated by the cosmological constant. Since we use an exact framework instead of first-order 
linear perturbations, as is done with FLRW models, our work is not limited to the linear regime or only first-order 
terms in the density contrast. We derive differential equations for the density contrast in both classes (converting to 
the growth rate variable G in class II), as well as expressions for the shear, expansion, and tidal gravitational field in 
terms of the density contrast and practical measurable cosmological parameters. In this form, the equations provide 
a more straightforward connection between theory and observation than do the metric functions on their own. We 
treat the class I formalism in much the same way as in [20j by introducing a quasilocal average density [36l |40| - |42T | in 
order to define a density contrast. To compare Szekeres growth to observations, we use data for the growth factor / 
from redshift space distortions and the Lyman-a forest to obtain best fit parameters for class II models and compare 
them to those of the concordance ACDM cosmology. We also explore the applicability of the usual growth index 7 to 
parametrize the growth factor in Szekeres. 

The outline of the paper is as follows. We first present in Sec. [TT| the GW formulation of the Szekeres solution. 
We then derive the exact growth equations for class I and class II models including a cosmological constant in Sec. 
IIIII Results and discussion from integrating the equations are given in Sec. IIV1 and in Sec. |V]we examine the early- 
and late-time behaviors of scalar quantities that appear in the Raychaudhuri equation and can be identified as the 
cause of the strengthened Szekeres growth. We derive the equation for the growth factor in class II models in Sec. IVII 
and compare it to cosmological data along with the standard ACDM cosmology. We conclude in Sec. IVIII Units are 
chosen so that 8irG = c = 1 throughout the paper. 



II. THE SZEKERES COSMOLOGICAL MODELS 



We begin with an introduction of the Szekeres exact solution to Einstein's field equations in the form due to Goode 
and Wainwright (GW) [HH 

and with a cosmological constant [35[ . The GW form is a reformulation of the original 
solution discovered by Szekeres [l8l H9|. A third formulation, discussed briefly in Appendix B, is also common in the 
literature and has a form similar to that of the LTB metric. It has been used in, for example, [27l I29T l30l . I32M34I l39l. |43| . 
The GW formulation is well suited for studying structure growth, since, as was remarked in |23f . one can consider 
the models as nonlinear exact perturbations of some averaged background. This background can be associated with 
a corresponding unperturbed FLRW model. The evolution equations in the Szekeres models can then be compared 
to linear perturbations of this associated FLRW model. It is worth noting that this association between the Szekeres 
background and an FLRW model is not the same as the case where a given inhomogeneous model becomes an FLRW 
model when some of the metric parameters (or certain coordinates) are taken to some large limit. One has to make 
further explicit definitions of some of the metric functions in order to implement such limits if needed. The full GW 
formulation (without a cosmological constant) can be found in the original papers [22|, [23j , and we therefore present 
only a brief introduction in order to be self-contained and to set the notation. The GW form of the Szekeres metric is 

ds 2 = -dt 2 + S 2 [e 2u (dx 2 + df) + H 2 W 2 dr 2 ] , (I) 

where it is assumed that the functions S, H, and W are all positive. The coordinates of the metric are comoving and 
synchronous, so the cosmic dust fluid has four- velocity components u a = Sq . 

We note that the coordinates x and y are not the "Cartesian" coordinates x and y. We have added a tilde on the 
GW coordinates for the Szekeres models in order to make this distinction clear to the reader. The GW coordinates x 
and y are the result of stereographic projections, and their transformation relations to known coordinate systems are 
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given in Appendix B, see also [43j. Next, for clarity we labeled the third GW spacelike coordinate as r instead of the 
z that was used in previous papers |2(j, H2, [23| . 

From the outset, we will refer to the scaling function of the spatial part of the metric as a instead of S to facilitate 
comparisons with ACDM. The metric function H depends on all four coordinates and can be written as the difference 
of two functions, only one of which carries the time dependence, 

H(t,r,x,y)=A{r,x,y)-F(t,r). (2) 

The function F satisfies the second-order linear differential equation 



P + 2±F-^-F = 0, (3) 

where ' = d/dt and M(r) is an arbitrary but sufficiently smooth function. This equation can be derived from either the 
field equations or from the Raychaudhuri equation for irrotational dust p3l l44j . There are two linearly independent 
solutions to Eq. §5§ denoted / + and /_ (the so-called growing and decaying modes, respectively), and so F can be 
written generally as 

F = /3+/+ +/?_/_, (4) 

where /3± are functions of r. 

The models divide naturally into two classes according to the metric function dependencies, and the exact forms of 
the functions A, v, f±, /3±, etc. for each class are given in Appendix A. In class I, which is the more general of the 
two, a = a(t,r), f± = f±(t,r), M — M(r), and W — W(r). In class II, these functions lose their r dependence so 
that a = a(t), f± — f±(t), M = const, and W = 1. 

The scale function obeys a generalized Friedmann equation 

a 2 _2M A k 

^~^ + 3"~^' ^' 

where k — 0, ±1. Equations Q and ([3]) govern the time evolution of the models and apply generally to both classes. 
By Eq. ([3]), we see that the class I a(t, r) satisfies the usual Friedmann equation of FLRW for every value of r so that 
the surfaces of constant r evolve independently in time. 
The matter density in both classes is given by 

6MA 6M ( F\ 

In class I, we can identify an exact density contrast S = F/H, which measures deviations from a background density. 
However, the interpretation of this 5 is different from that of the usual FLRW S — (p — p)/p, where p is the average 
density of the space. S at some event (to, tq, Jo, ya) compares the density there to the average density inside the surface 
defined by to and ro, instead of to some overall background or limiting value of p where the solution becomes FLRW. 
Therefore, in order to facilitate the comparison between the growth in Szekeres class I and the perturbed FLRW, we 
use a Szekeres model as specified in Appendix C, representing a large-scale cosmic structure surrounded by an FLRW 
model. On the other hand, in class II, a and M have no r dependence, so it is natural to define a background density 
by p(t) — 6M/a 3 (t). By keeping 5 = F/H (no hat to distinguish it from class I), the density contrast then describes 
deviations from some smooth underlying background like it does in FLRW. But while S is written in the same way as 
in FLRW, a spatial profile is still needed in the Szekeres class II in order to facilitate the comparison to the FLRW 
plus perturbations scheme. This can be done for class II using the metric functions f3±(r). Indeed, one recalls that 
both classes of Szekeres solution become FLRW when the functions f3±(r) are zero. This is in fact the necessary and 
sufficient condition 23]. It is in this case, that F and 5 become identically zero, and the density assumes a similar 
form, p = 6M/a 3 , to that of FLRW. Since class I has been claimed to be more relevant than class II for astrophysical 
applications [43| . we chose to implement the profile modeling in class I only and use the function M(r), while for class 
II, we integrate the growth over time for a fixed value of r with no further modeling of the functions j3±(r). Finally, 
we note that the density is constrained to have the same sign as M, and so to have a physically reasonable solution 
we restrict our investigation in this work to models with positive M. 



III. COSMOLOGICAL EVOLUTION OF LARGE-SCALE GROWTH USING THE SZEKERES MODELS 



In class II, the density contrast offers a straightforward comparison between exact Szekeres models and the linearly 
perturbed FLRW equation in 6. However, we first address large-scale structure growth in class I via 5, as it is the 
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more general of the two classes and requires more care. The treatment here is similar to our previous work [20| and 
employs, for class I, the quasilocal variables introduced and used in (36l. l40l - |42j . The main difference between this 
and our previous work is that whereas before we considered only the matter-dominated era, we now derive growth 
equations that include a cosmological constant, and so our analyses apply from the time of last scattering up to the 
A-dominat ed p resent and future evolutions. Finally, we recall that the Szekeres class I flat models have no growing 
modes [13, Elf , so we left them out. 



A. Growth equations in curved Szekeres class I models with a cosmological constant 



In order to treat the two Szekeres classes similarly, we make some definitions for class I so that the matter density 
may be split into a background and a density contrast, as mentioned at the end of Sec. [TTJ We introduce the quasilocal 
average density [H, l40l - |42j 



Pq(t,r) 



fy fx Ir Fptf, r, x, y)V^h dr dx dy 



(p)qV[r](t) 



(7) 



over the domain T>[r], which is bounded by an r = const surface and is a subset of the three-dimensional hypersurfaces 
of constant t. The function J- is db weighting factor, and the projection tensor {h^y — ~\~ u^Ui, ) has determinant 
h of its three-dimensional part. The density contrast 6 can then be defined as usual but with p q {t, r) serving as the 
background average, 



p{t,r,x,y) - p q (t,r) 



Pg{t,r) 

As discussed in [3(1 HoMHI , p q is a coordinate invariant quantity that can be shown to be given by 

6M(r) 



P q (t,r) 



a 3 (t,r)' 



By Eq. ([8|), we then have 
which, upon comparing to Eq. 



p(t,r,x,y) = p 9 (t,r)[l + 6(t,r,x,y)], 
allows us to make the identification 

F(t,r) 



H(t,r,x,y)' 



(8) 



(9) 



(10) 



(11) 



Since p q is a coordinate invariant quantity, so too must 8 be. As was done in [20], the time differential equation for 
F [Eq. ([3])] can be recast in terms of 6 to arrive at the class I evolution equation for the density contrast, 



a(t,r) a 3 {t 7 r) \ + § 



a 3 (i,r) 



0. 



(12) 



The form of the above equation is suggestive of linear density perturbations in FLRW. To clarify the connection, we 
use Eq. (J5J) to define generalized density parameters analogous to those of ACDM. Such a generalization is similar to 
what has been done in previous works for the LTB inhomogeneous models in order to compare them to observations 
(e.g., see [l3l - [T5j and references therein). Equation ([S]), with its explicit coordinate dependencies, in class I is 



« 2 (^) = m , = 2M(r) A _ k 
a 2 (t,r) ~ 1 ' ; a 3 (t,r) 3 a 2 (t,r)' 



(13) 



where we have defined H by analogy with the Hubble parameter in ACDM. We point out that H is not the same as 
the metric function H . Continuing by analogy, the t- and r-dependent Szekeres cosmological density parameters can 
then be defined as 



Q m (t,r) 



2M(r) 



n A (t,r) = - 



a 3 (t,r)M 2 (t,r) 
A 



3H 2 (t/, 



(14) 
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and 

so that Eq. (fl3| can be written as the sum of the three equal to 1, as is usual in ACDM. 
By Eq. © , we can now write Eq. (TT2"j) as 

5 + 2M(t,r) § - 4irG Pq (t,r) 6 S 2 - 4nGp q (t, r) S 2 = 0, (15) 

1 + S 



where we have temporarily restored the factor 8nG for clarity. Equation (|15|) . in its first three terms, bears a formal 
similarity to the linearly perturbed FLRW equation in the density contrast. However, the Szekeres 5 measures 
deviations of p from p q , the quasilocal average density, as opposed to a global average as in FLRW. The H and p q 
terms are also dependent on r, and Eq. (fT5"|) contains two nonlinear terms that the FLRW counterpart does not. 
Moreover, 5 here is an exact quantity arising directly from an inhomogeneous metric and as such is not constrained 
to be smaller than 1. 

In class I, the evolution equation for 5 in terms of time derivatives is most useful, but we can continue the formalism 
in a general way that will be more useful for class II, where a derivatives and density parameters are convenient. Since 
each surface of constant t and r evolves independently, we can consider fixing the value of r at some rg and letting 
a{t, ro) serve as the time parameter for that surface. Then in terms of a derivatives, Eq. (|12p becomes 

A + l>-^-^-3^ = 0) (16) 

a z a J a 3 a z 1 + 5 a J a z 
where ' = d/da at ro. 

An explicit expression a (parametric or otherwise) is unnecessary, since we can eliminate d and a in favor of Sl m 
and SI a, using Eq. (fl~3|) and its time derivative. Multiplying Eq. (fl"3| by a 2 and differentiating with respect to t gives 

MA 

a = -- + -a. (17) 



With this and the definitions (O, we obtain 

a 2 4 + 2S1 A - SI,, 



■,2 



and 



a 2a 
3Af 3 Sl m 



(18) 



2 a 



( 19 ) 



Substituting these into Eq. P^|) . we find 

-„ U + 2 S1 A - Sl m \ _ 3 _ _^y 2 _ 3 ilrn^ = Q 

V 2a / 2 a 2 i + S 2 a 2 

Finally, if we assume the cosmological parameters take on their standard meaning at every r, we can express them 
in terms of their values evaluated today (denoted by a superscript naught) and a(t, r), 

O (t r) = n< L( r ) (2i) 

Hm{h r) QP m (r) + n° (r)(a/oo)3 + fig(r)(a/ao) 1 j 

and 

S1 A (t r) — 0^(0(^)3 

where SI? (r) = 1 — n^(r) — Sl A (r) and ao = a(to, r) is the scale factor today. These relationships can then be evaluated 
at ro and substituted into Eq. (|2U|) . 

In class I, there is an alternative method of calculating 5 when a model is completely specified — that is, when all 
the arbitrary functions of r are known (see Appendix A) . Instead of numerically solving the second-order differential 
equation in (|12[) . one can compute 6 via the explicit functional form 

, M, r a ir \ ( a, 



where a comma denotes partial differentiation. 
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B. Growth equations in flat and curved Szekeres class II models with a cosmological constant 

As discussed in Sec. [TTI above, the class II density contrast 5 = F/H=(p — p)/p has a similar form as in FLRW, 
where p(t) = 6M/a 3 (t) is the overall background average density. We also understand pit) as the associated FLRW 
model's density that underlies the exact Szekeres perturbations. This is possible since a and M no longer depend on 
r. Equation (j2"0)l applies equally well to class II but with 5 replaced by 8, a = a(t), and Hi = fi^(i) for i = {m, A, k}. 
Now the density parameters in terms of their values today are 

fi° 

and 

= o ? „ + ntW - (25) 

where fi° = 1 — fi^ — fi^v and we have set ao = a(to) = 1. In class II we can therefore directly compare observations 
with the cosmological density parameters between Szekeres and ACDM. 

As in our previous work [201 ] . we find it useful to rewrite the evolution equation for d in terms of the growth rate 
variable G = (S/5 )/a. Noting the relations 6' = G + aG' and 6" = 2G' + aG" (for the choice S = 1), Eq. (J2DJ for 
class II then becomes 

G" + f 4 + fi A %i) + (2 + n A - 2 n m )% - - {G + aG ' f - §fi ro — = 0, (26) 
\ I J a a z a 1 + «Cr 2 a 

and one can easily check that this equation reduces to the forms obtained in 20] for the special cases (i) fi A = and 

(ii) n m = i, n A = n fc = o. 



IV. INTEGRATION OF GROWTH HISTORY FOR SZEKERES MODELS WITH A COSMOLOGICAL 

CONSTANT 

The Szekeres growth equations derived in Sec. IIIII are integrated numerically for 6 (class I) and G (class II) using a 
standard fourth-order Runge-Kutta routine with adaptive step size [45] . In class I we discuss results for a simple but 
inhomogeneous and nonsymmetric model representing a central void with neighboring supercluster, which is similar 
to Model 1 employed by [2{|- Our model is designed to match to a nearly homogeneous FLRW background model 
described with standard ACDM parameters (fi°„ = 0.27, = 0.73, and H = 72 km s _1 Mpc -1 ) at a radius of 50 
Mpc and greater. We provide the algorithm for building such a model, specify the necessary arbitrary functions, and 
describe the resulting density structure in Appendix C. 

For class II, we compare the growth rate in Szekeres models with various combinations of cosmological parameters 
to the ACDM model with linear perturbations. The discussions for class II can be applied as well to class I if we 
restrict ourselves to a constant r-coordinate surface, where the density parameters take on the values specified for the 
class II model. However, as seen by Eq. (|A7j) of Appendix A, there are no growing modes for the spatially flat class 
I case, because k = forces /3+ to vanish. Since we are interested in structure growth, we therefore do not consider 
flat models to apply to class I. 



A. Growth history for flat Szekeres class II models 

We find that for flat class II Szekeres and ACDM models with comparable values and 0° , the Szekeres models 
exhibit significantly stronger growth than their linearly perturbed ACDM counterparts. This difference is as expected, 
since the Szekeres growth equation contains exact nonlinear contributions due to inhomogeneities as well as shear and 
a gravitational tidal field that the standard concordance model does not. 

The left panel of Fig. [T] shows the growth behavior over a range of f2° n in flat models. For comparison, the linearly 
perturbed ACDM growth for = fig + fi^ m = 0.27 and Q° A = 0.73 is plotted as well. The Szekeres case with 
fi^ = 1.0 and ft A = (i.e., with an Einstein-de Sitter associated background) experiences the strongest growth with 
no sign of suppression within our cosmic history. For smaller values of fi^, however, we begin to see suppression, 
which is consistent with the fact that the cosmological constant makes the expansion of the universe accelerate and 
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FIG. 1: LEFT: Growth rate of large-scale structure in flat Szekeres class II models with a cosmological constant. A matter-only 
universe experiences the strongest growth, and curves are seen to experience significant suppression by A within our cosmic 
history for fij^ values of approximately 0.1 and smaller. The curve with fi^ = fig = 0.04 and fii = 0.96 ( green, solid) could 
provide a consistent growth history and is a good flat model candidate for comparison with ACDM (black, dashed), which is 
plotted with the standard values of fi^ = 0.27 and fi° = 0.73. RIGHT: A sample of flat class II models guided by the results of 
the left panel that could mimic ACDM (dashed) with respect to the growth of large-scale structure today. The Szekeres curves 
emerge at early times almost identically but begin to separate at a scale factor of about 0.15. The fig = 0.05 model is able to 
grow the most structure before suppression by the cosmological constant takes over. Despite differing times for the onset of 
suppression, the three curves experience approximately the same rate of suppression at later times, indicating the effect of a 
dominant cosmological constant. 



gives large-scale structure less opportunity to grow. In each case, the cosmological constant contribution ultimately 
wins out, turning the curve over, but only for VL m ss 0.1 and smaller do we see the effect clearly in our cosmic history. 

It is interesting that if we take into consideration the bounds 0.039 < fig < 0.049 from big bang nucleosynthesis 
(BBN) (4fj| . we can interpret the small matter contribution in Szekeres as being due solely to baryons (i.e., not dark 
matter). With the BBN limits in mind, we find that a flat Szekeres model with = 0° w 0.04 produces a growth 
curve that roughly mimics that of linearly perturbed ACDM today. The Szekeres growth rate is larger early on than 
ACDM but still undergoes suppression to make 5 /a of the two comparable at present. This is significant because 
considering the effect of nonlinear inhomogeneities seems to strengthen the growth of structure in a way that requires 
little or perhaps no dark matter. We discuss this point further in Sec. IVI1 where we compare the Szekeres growth to 
current observational data and obtain best fits for the parameters O^, , and Sl^ . 

In the right panel of Fig. [TJ we see indeed that flat Szekeres models with values of f2^ = 0° near 0.04 have growth 
rates comparable to ACDM today, with the spread arising from shifting the balance between and In each 
curve here, the growth rate emerges at early times almost identically with the others up to a ~ 0.15, but the curves 
representing larger values experience suppression later. As expected f47| . when A dominates at late times, the 
growth rate becomes essentially linear with constant (negative) slope among the three cases presented. 

B. Growth history for curved Szekeres models 

While the global spatial curvature of the universe has been well constrained to be negligibly small using the ACDM 
model [46} this is not the case when using inhomogeneous cosmological models [48j ]. As we will see in Sec. IVI1 the 
zero curvature result using FLRW may be biased due to assumptions that do not hold in the more general Szekeres 
models. We compare there class II models to growth data and find a significant curvature component, which is absent 
in the ACDM model. 



1. Class I 

The evolution of the density contrast for our class I model (described in Appendix C) is shown in Fig. [3] Though 
the model matches (at some given r) to a nearly homogenous ACDM background with zero curvature (f2^, = 0.27, 
= 0.73), the Szekeres region itself has a nonzero curvature (k = —1). The left panel shows 5 as a function of r 



r [Mpc] 



r [Mpcl 



FIG. 2: LEFT: Evolution of the density contrast of the class I model along a direction through the overdense region that 
neighbors the central void. The color gradient represents 8 at different times, where the present time is darkest and curves get 
lighter toward the past. RIGHT: Evolution of the density contrast along a direction orthogonal to the overdensity. While in 
the left panel S peaks as a function of r, in the right panel it is a monotonically decreasing function of r that never exceeds 1. 
In both panels, 5 vanishes beyond 50 Mpc, which we expect due to our model matching to an almost FLRW background at 
that radius. The curves are evenly spaced in time, and the density contrast in both directions is seen to decrease toward the 
past, indicating a smoothing out of the region. 



in a direction passing through the overdensity adjacent to the central void. There is a prominent maximum today at 
r f=a 11 Mpc, and this peak diminishes toward the past as well as moves toward r = 0. The density contrast profile in 
this direction achieves the largest values compared with other directions, which is consistent with the expectation that 
the density deviates most from the background through regions of larger structure. However, we must keep in mind 
that the background density at a given r is a quasi-local average density p q defined as the average within the volume 
bounded by the 2-surface of constant r and t. The maximum of 5 therefore need not coincide with the maximum 
of p for some value of t (see Fig. [5]). The right panel plots 5 along a direction orthogonal to the overdensity. The 
density contrast here is a monotonically decreasing function of r for all times shown, and the density deviation from 
the background is seen to lessen toward the past. It is clear that less structure lies along this direction compared 
to that of the left panel. That our model matches to an almost FLRW background (as described at the beginning 
of Sec. IIVP at 50 Mpc is apparent in both plots, since 5 approaches zero there. We find here that the growth in a 
Szekeres class I model is up to several times stronger than that of the linearly perturbed ACDM model (S <C 1), a 
result consistent with previous works (2p| . |30] | . 

2. Class II 

The inclusion of curvature into the growth for class II provides another degree of freedom that allows for a variety 
of different growth histories. To see clearly the effect of all three components (matter, curvature, and A) on S/a, we 
fix one component at a time and vary the other two. For example, in Fig. [3] we first set Q® n — 0.1 and vary and 
f2°. For the same matter content, closed models (fi^ < 0) experience stronger overall growth than open (O^ > 0) 
models, with the flat case falling in between. At early times, curvature seems to be the determining factor in the 
growth rate, since it is the curve with 17° = —0.02 (and also the largest value of ) that rises quickest early on. But 
as expected, later the cosmological constant effects its characteristic suppression so that the evolution of each curve 
is approximately the same after it reaches its maximum. 

Next, the effects of fixing at 0.9 and varying and are shown in the right panel of Fig. [3] This plot 
spans the same curvature range as the one to the left, yet here G exhibits more pronounced variation. In this case, 
closed models correspond to having more matter, while open models to having less. The trend in the curves therefore 
matches the expectation that more matter in a closed universe will have a stronger growth rate of structure than in 
the opposite case. Again, A suppresses the growth uniformly after the curves turn over, but now the maxima occur at 
very different times. We also note that the presence of nonzero f2° causes the curves to diverge from each other quite 
early — even earlier than a ~ 0.1. Flat models for the parameters in Fig. [TJ by contrast, evolve almost identically up 
to scale factors of about 0.15. 

As we will see in Sec. I VII we can find parameters such that the Szekeres growth history exactly matches that of 
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FIG. 3: LEFT: Class II Szekeres growth histories for 0°, = 0.1 and varying Q.° A and fl° k . Closed models (fi£ < 0) experience 
stronger growth than open models (f2^ > 0) overall. The flat case lies in between the two, and all models have stronger growth 
compared to standard ACDM (not shown; see Fig. [1} ■ All the Szekeres curves undergo the same characteristic suppression by A 
after attaining their maximum growth rates (which happens later the smaller is) so that the slopes becomes approximately 
the same across the different curves. RIGHT: Szekeres growth histories for Q% = 0.9 and varying and SI®. Closed models, 
which have larger matter density, again experience stronger growth than open models with the flat case lying in between. The 
spread in the curves is slightly more pronounced here than for fixed 0„ over the same range of fijj values, indicating the effect 
that small changes in the matter content can have. All models have stronger growth at all times than ACDM (not shown; see 
Fig. [T]) and ultimately get suppressed by A before a reaches 1. 



ACDM. The parameters that accomplish this are ilP rn = 0.11, S1° A = 0.71, and 0£ = 0.18. In Fig. [TJ its G curve 
would be indistinguishable from ACDM (black, dashed). In this case, if we follow the predictions of BBN, the need 
for dark matter cannot be fully eliminated, but the amount needed to achieve the same growth as ACDM is reduced 
by a factor of about 3. Also, here the dark energy density (or cosmological constant) contribution is nearly the same 
as for ACDM. The Szekeres ilP k is significantly higher than that of ACDM (which is consistent with zero). The Si's 
in the Szekeres models do represent fractions of the total energy content of the universe for different constituents, as 
defined in the Friedmann equations, which are formally identical in both the Szekeres and ACDM models. However, 
the energy densities do not have strictly the same connections to physical quantities in both cases. For example, the 
Hubble parameter in FLRW is proportional to the expansion scalar O, but in Szekeres there is also a dependence on 
the metric function H (see Sec. IV]) , 

Based on the results of integrating the growth equations, we find that class II Szekeres models do have the ability 
to mimic the growth rate of large-scale structure in ACDM today for appropriately chosen values of the cosmological 
parameters. In particular, the parameters il® n = 0° « 0.04 and f)9 ~ 0-95 (so that the universe has a small 
negative curvature of order 10 -2 ) do a fair job qualitatively, and interpreting the matter content as being solely due 
to baryons supports the idea that Szekeres models can provide a consistent growth history with little or no dark 
matter. Alternatively, the parameters = 0.11, — 0.71, and = 0.18 exactly reproduce the ACDM growth 
history, but seem to require some dark matter and large 0^. We present a quantitative argument for smaller $7^ in 
Szekeres than ACDM based on fits to the data for the growth factor in Sec. I VII 

Finally, we also observe that the same trends in the growth rate hold in general for FLRW models as they do for 
their Szekeres counterparts. In other words, closed universes with larger matter contributions experience stronger 
growth overall than open universes with less matter. Flat universes lie in between. The effect of a cosmological 
constant in ACDM is also to suppress the growth rate equally at late times when it becomes dominant, regardless of 
the growth rate at earlier times, a well-known and expected result, e.g. 



V. COSMOLOGICAL EVOLUTION OF SHEAR AND TIDAL FIELD SCALARS IN SZEKERES 
MODELS INCLUDING CURVATURE AND COSMOLOGICAL CONSTANT 



We now consider the time evolution of physical scalar quantities in the Szekeres models that influence the gravita- 
tional clustering and structure growth. Our discussion here is brief, since much of the work was carried out in our 
earlier paper [20j , and the inclusion of A in the expressions for the shear scalar and other quantities is straightforward. 
In class I, due to the r dependence of a, it is natural to express these quantities in terms of 5, while in class II we 
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FIG. 4: LEFT: Evolution of |<x/|0 in the class I model along a direction through the overdense region that neighbors the central 
void. The color gradient represents <t/|0 at different times where the present time is darkest and curves get lighter toward the 
past. The maximum is achieved today at r w 20 Mpc, and the shear-to-expansion ratio diminishes continuously toward the 
past with the peak moving toward r = 0. RIGHT: Evolution of |<r/|0 in the class I model along a direction orthogonal to the 
overdensity. The curves in this case first increase as we move backward in time before ultimately decreasing as in the left plot. 
The curves are evenly spaced in time, and so | a - / 1 © is changing at a slower rate today than it was in the past. We expect the 
shear to tend toward zero in the past as the model smooths out, and looking through the structure in both directions reveals 
that this is the case. 



convert to G as in our previous work. 



A. Shear and expansion 



The first quantity of interest is the squared shear scalar defined by 



V^ 1 fl ) 



where the mixed tensor components are given by |49j 

2<j x x — 2cr y y = —o r 

We can write then for class I, 



2H 
3H' 



(27) 



(28) 



and for class II, after some algebra, 



2 S 2 



3(l + 5) 2 



G + aG' 
1 + aG 



(29) 



(30) 



where Ho is the Hubble parameter today and do = 1- As we noted in previous work [20j , nonzero shear augments the 
gravitational attraction, resulting in stronger collapse than would be produced by the energy density alone. 

We next give the expression for the expansion scalar in order to examine early- and late-time divergences of 
| <7 1/0. The expansion scalar is denned as 



= 



a H 
a H 



(31) 



which we find for class I is simply 
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FIG. 5: LEFT: Shear scalar over the expansion for class II Szekeres models with no cosmological constant. We see |c|/G going 
to zero at early times and potentially diverging as a becomes large depending on the amount of matter and curvature. The flat 
matter-only case diverges most strongly, and the curves grow without bound down to f2„ ~ 0.4, below which we see the curves 
turn over and remain finite in the future. RIGHT: Shear scalar over the expansion for class II models with a cosmological 
constant and fixed Q.^ = 0.1. \cr\/Q in all these cases goes to zero for small a as in the models of the left panel. The curves 
reach a maximum between 0.4 < a < 0.8 and then tend to zero again as a becomes large with no divergences. In models where 
fl% > 0.9 (positively curved), the ratio |cr|/© attains a higher maximum at later a than those where < 0.9 (negatively 
curved), and the flat case falls in between. 



and for class II can be written 

e = B °(l + ^)^ + n > 2 + a i- < 33 > 

We plot |ct|/6 for the class I model described in Appendix C along a direction passing through the overdensity 
(left) and one orthogonal to it (right) in Fig. 2) The left plot has its maximum curve today with a peak at r « 20 
Mpc, and the curves there decrease with peaks shifting toward smaller r as we move toward the past. In the right 
plot, where we are looking through less structure, the peak in |er|/0 today also occurs at about 20 Mpc, but toward 
the past now the curves first rise before ultimately decreasing. In both cases, the ratio of shear to expansion vanishes 
at r > 50 Mpc, which is expected due to the matching: all FLRW models have zero shear, and the one we match to 
has nonzero (positive) expansion. Also, since the curves are evenly spaced in time, we see that |cr|/0 was changing 
faster in the past than it is today for both directions, and nowhere does it diverge for the times and ranges plotted. 

The left panel of Fig. [5] shows |oj/0 for various class II models without a cosmological constant. These and the 
following class II figures have been extended beyond a = 1 to see the trends in the curves more clearly. The ratio 
of shear to expansion goes to zero as a approaches zero, indicating an isotropic initial singularity, and it potentially 
grows without bound in the future depending on the values of the matter and curvature density parameters. We 
can also see this by taking limits of the expression for |tr|/0. The matter-only case grows most rapidly and shows 
no signs of turning over down to values of approximately 0.4. Below this value, where the models are strongly 
negatively curved, the shear and expansion evolve in proportion to each other, and the ratio does not diverge. The 
right panel shows |oj/0 for class II models with a cosmological constant and fixed at 0.1. Again, all curves shown 
go to zero for small a but now also tend to zero as a becomes large, attaining a maximum between scale factor values 
0.4 < a < 0.8. Curves for models with f2j^ > 0.9 (positively curved) reach higher maxima and at a later a value than 
those where < 0.9 (negatively curved), and the flat case falls in between. 

We note that the explicit 51° dependence cancels out in the expression for |er|/0 in class II, which makes it identical 
to that previously found for Szekeres models without a cosmological constant. However, is bound up in the 
equation for G, so nonzero fi^v still affects the growth history and therefore |er|/0 as well. This same observation 
applies to the density p and the tidal gravitational field scalar E that we examine next. 
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FIG. 6: LEFT: Evolution of the density (normalized by the almost FLRW background density pb) of the class I model along 
a direction through the overdense region that neighbors the central void. The color gradient represents p/pb at different times 
where the present time is darkest and curves get lighter toward the past. This is the region with maximum density in the 
model, but it smooths out (i.e., becomes more homogeneous) significantly in the past, as expected. RIGHT: Evolution of the 
normalized density along a direction orthogonal to the supercluster. The maximum overdensity attained here is about half 
that of the left panel, while the underdense regions in both have approximately the same profiles at corresponding times. The 
model also smooths out in this direction toward the past. In both panels, the curves are evenly spaced in time, from which we 
can infer that though the rate of growth of the overdense region is undiminished up to the present, the rate of growth in the 
regions with less structure seems to decrease with time. The central void's volume continues to grow in both directions over 
time in order to compensate the growth of the overdensity. 



B. Density and tidal gravitational field 

As seen in Sec. IIII1 the equation for the density in class I is given by 

P=— (1 + 5), 34 

or 

which can be written in terms of G for class II as 

3f7° 

p= ±^±B 2 (l + aG), (35) 

where we have used the definition of Q m - 

In Fig. [S]we plot the density, normalized by the almost FLRW background density pb, of our void plus supercluster 
class I model as a function of r, again looking through (left panel) and orthogonal to (right panel) the overdensity. As 
expected, the region with the supercluster has the highest density today, and in both panels the deviation from the 
background becomes smaller toward the past. If viewed in terms of an associated density contrast, where we consider 
density differences with respect to the limiting FLRW value, then we see the quantity p/pb — l go to zero as we move 
backward in time. That is, the overdensity and the void both get smaller and smooth out, leaving just a void with 
ever-shrinking volume. It is interesting that the rate of growth of the overdense region is undiminished up to the 
present, while the rate of growth of the regions with less structure seem to decrease with increasing time. However, 
the central void's volume continues to grow in both directions over time in order to compensate for the growth of the 
overdensity. 

The left panel of Fig. [7]shows /o/Hq for class II models without a cosmological constant. As can be understood from 
Eq. Q35p , all models there diverge as a approaches the initial singularity, but late-time divergences within our cosmic 
history depend on the values of the matter and curvature parameters. The matter-only case diverges at a « 0.5, and 
models with less matter (larger f2°) also have densities that diverge at later a with smaller tt® n resulting in a later 
onset of divergence. The right panel shows the evolution of density in Szekeres models that include a cosmological 
constant and have set to 0.1. The models here also all diverge at the initial singularity, but they remain finite 
up to today and tend to zero in the future, decreasing monotonically. Models with 17° > 0-9 (positively curved) have 
density curves lying above those with < 0.9 (negatively curved), and the flat case falls in between. 

We now turn to the tidal gravitational field, and for that we consider the electric part of the Weyl conformal 
curvature tensor. Its components are defined as 



(36) 
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FIG. 7: LEFT: Evolution of density in class II Szekeres models with no cosmological constant. As expected, all models here 
diverge as a approaches the initial singularity, but late-time divergences within our cosmic history happen at times that depend 
on the combination of matter and curvature parameters. For example, the ft^ = 1 case diverges at a ~ 1, followed by the 
model with Sl^ = 0.8, and models with less matter experience density divergences in the future at later a values. RIGHT: 
Evolution of density in Szekeres models with a cosmological constant and fixed Q.^ =0.1. Again, all models diverge at the 
initial singularity, but they remain finite up to today and tend to zero as a becomes large. Positively curved models (Q% > 0.9) 
have density curves lying above the negatively curved models < 0.9) with the flat case in between. The addition of a 
cosmological constant avoids the late-time divergences in the density. 



where the Weyl tensor in four dimensions has components 

1 R 

C^vufi = Rfivafi + -^(g^pRua + 9vaRp,p ~ 9^ a Rvfi ~ 9vf$Rva) + -^{9^a9uP ~ 9jifi9va)- (37) 

The magnetic part of the Weyl tensor is zero in the Szekeres solution, and the models therefore fall into the 
category of "silent" cosmological models. The shear and the electric part of the Weyl tensor are related via the 
propagation equation for the shear 

+ x<T X „ + |e<7" „ - „ + it" u v )a 2 = -W „, (38) 



which reveals how E llVl interpreted as the tidal gravitational field [17|,l50|, causes shearing in the fluid flow. 
We are interested in how the tidal field scalar, the magnitude of which is given by 

E 2 = E%E\, (39) 

and whose mixed components are [5l| 

2^. = 2^ # = -^ r = |f| + 2^§) J (40) 



compares to the energy density at early and late cosmic times. The class I expression has the simple form 

E 2 =&(^5]\ (41) 



a° 

and in class II, in terms of cosmological parameters evaluated today, a small calculation gives 

B 2 = |(^G) 2 . (42) 

Again, in the expressions for p and E 2 , the dependence on the cosmological constant comes via a and 17° , which 
is contained within G. Figure [S] shows the ratio \E\/p along the two directions considered previously for our class 
I model described in Appendix C. The left panel is a radial profile through the supercluster, and the right panel is 
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FIG. 8: LEFT: Evolution of \E\/p for the class I model along a direction that passes through the supercluster neighboring the 
central void. The color gradient represents \E\/p at different times where the present time is darkest and curves get lighter 
toward the past. It evolves from a localized peak around r ~ 1 Mpc in the past to a larger and significantly broader maximum 
around r w 11 Mpc today. RIGHT: Evolution of \E\/p along a direction orthogonal to the supercluster. Here, the tidal 
field-to-density ratio is a monotonically decreasing function of r that diminishes continuously toward the past. In both plots, 
\E\/p approaches zero at r = 50 Mpc, which we expect due to matching to the almost FLRW background there and that the 
Weyl tensor is zero in FLRW models. In addition, as adjacent curves are spaced equally in time, the tidal field was changing 
more quickly with respect to the density in the past than it is today. 




FIG. 9: LEFT: The ratio of the tidal gravitational field scalar to the density in class II Szekeres models with no cosmological 
constant. The curves all grow monotonically from zero at the initial singularity and remain finite within our cosmic history. 
The larger the matter content, the larger the ratio of \E\ to p, and the less suppression the curve experiences at later times. 
RIGHT: The ratio of the tidal gravitational field scalar to the density in class II models with a cosmological constant and 
fixed = 0.1. As in the case without Q,a, all the curves grow monotonically from zero, but here they all also evolve almost 
identically after a ~ 1. As a becomes larger, |_E| and p for a given set of parameters ultimately grow at the same rate so that 
their ratio approaches a constant. The larger curves have larger \E\/p than those with smaller Q% for the same matter 
content. Again, the addition of a cosmological constant removes the late-time divergences seen to the left. 



plotted for an orthogonal direction. In both cases the tidal field-to-density ratio increases continuously up to the 
present, but in the direction of the overdensity, the curves evolve toward wide peak around r« 11 Mpc. In the other 
direction, the curves are monotonic decreasing functions of r. The matching to the almost FLRW background at 
r = 50 is again clear in both panels, which we expect since the Weyl tensor is zero in FLRW models. Also, since the 
curves are spaced evenly in time, we see that \E\ was changing more rapidly compared to p in the past than it is at 
present. These observations are reminiscent of the behavior of 6 itself, which is not surprising, because \E\/p can be 
written in terms of only 5 as \E\/p — <5/[\/6(l + §)]. 

We plot \E\/p in Fig. [5] and compare early- and late-time divergences in class II models with A to those without. 
On the left are models without a cosmological constant, and we see that all curves grow monotonically from zero at 
a = and remain finite within our cosmic history. As with |cr|/0 and p/Hq, the matter-only case is on top, and larger 
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TABLE I: Summary of observational data used to constrain the cosmological parameters in the Szekeres class II and ACDM 
models through the growth factor. Data is taken from the compilation of [52J and [53(], with original references for each value 
listed in the right column. The first nine values are derived from measurements of redshift space distortions, while the last seven 
are from Lyman-a measurements. Where necessary, we have given the mean redshift of the / measurements when originally 
reported for a range of redshifts. 

matter content corresponds to a larger ratio of \E\ to p. On the right we plot Szekeres models with a cosmological 
constant and fixed = 0.1. As in the case without all the curves grow monotonically from zero, but here they 
all evolve almost identically after a 1, indicating that \E\ and p grow proportionally at larger a. Again, similar to 
| cr |/0 and p/Hq, larger curves overall have larger \E\/p than those with smaller for the same matter content, 
with the flat (f2° = 0.9) case falling in between. The important point is that for all the models plotted here, the 
addition of a cosmological constant removes any late-time divergences that were present without it. 

Furthermore, since the cosmological constant will continue to dominate as time increases, one expects a future 
evolution free from singularities for the models of interest. The Szekeres solutions can have shell-crossing singularities 
when the metric function H vanishes. This is clear from Eq. © for the matter density. We have verified that for 
the models of interest explored above with < 0.1, the growth rate and the density do not become singular in the 
future. G approaches zero as a becomes large, and since the scale factor grows without bound in these models, the 
density goes to zero along with G. 



In the preceding sections, we have demonstrated that the Szekeres class II models are capable of reproducing the 
required features of the growth history of large-scale structure using a nonzero cosmological constant A to mimic the 
ACDM concordance model. This includes both a late-time suppression of growth due to A and significantly stronger 
growth at early times than is possible with linearly perturbed ACDM using only a baryonic matter component 
0° = 0.04 consistent with constraints from Big Bang nucleosynthesis (BBN). While we were able to explore an initial 
parameter space for the Szekeres models by qualitative comparison to ACDM, we will now compare the models directly 
to currently available growth data. 

The following discussion only applies to class II models. We do not attempt to compare class I models to available 
growth data, as this requires a full consideration of the cosmological parameters as well as the available data as 
functions of r. There is not enough growth data points at the present time in order to attempt such r-modelization 
and we leave this for future work while specializing the current comparison to class II. 

We will use measurements of the growth factor 



which has been constrained over a wide range of redshifts. The values of / that we compare to are listed in Table 
U and come from measurements of both redshift space distortions and the Lyman-a forest. The growth factor is 
derived from either the redshift distortion parameter /? = //& (where b is the deterministic galaxy bias) or from the 
amplitudes of power spectra from Lyman-a forest data. Certain assumptions based on ACDM are necessary to obtain 
this growth data for / and thus may impact a direct comparison of the data to the growth resulting from modeling 



VI. CLASS II COMPARISON TO CURRENT GROWTH DATA 



/ = dlni5/dlna, 



(43) 
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La? 


Without priors 

fiX fi£ (x 2 ) 


With priors 

n£ (x 2 ) 


Szekeres y 


0.12 0.59 0.29 (0.22) 
0.11 0.69 0.20 (0.39) 


0.05 0.98 -0.03 (0.62) 
0.05 0.98 -0.03 (0.63) 


ACDM ^° 
Yes 


0.29 0.56 0.15 (0.21) 
0.26 0.69 0.05 (0.39) 


0.27 0.73 0.00 (0.32) 
0.27 0.73 0.00 (0.43) 



TABLE II: Best fit parameters for Szekeres and ACDM models. We use the \ 2 minimization described in Eq. (|52[) in order to 
choose the best cosmological parameters that fit / to the observed growth factor as given in Table [T] We present results both 
including and omitting available Lyman-a data. We also present results with priors placed on the Szekeres model, such that 
0.039 < fifj < 0.049 from BBN (with fi° m = 0). For comparison, we include the \ 2 f° r each best fit along with the x 2 f° r 
standard ACDM parameters. Values are rounded to the nearest percent. The Szekeres model with no priors fits the growth 
data with approximately the same \ 2 as that of the ACDM model. 



the evolution of structure in a Szekeres universe. However, we explore here the formalism of such a comparison and 
provide some first insight into the Szekeres models' ability to fit real growth data. 



A. Growth factor and growth index in the Szekeres models 

In ACDM, the growth factor is written from the growth equation as 

/'+^2+^/ + / 2 -^ m = ! (44) 

where H here is the ACDM Hubble parameter; see for example [52l [63 - lf38j . 

In a similar way to previous work done for ACDM, we can derive the growth factor for the Szekeres models. From 
the definition of / in Eq. f4"5|) . we substitute 

5' = f- (45) 
a 

and 

S"={f -/ + /') 4 ( 46 ) 
into Eq. (|20|) . Simplifying the resulting expression gives 

/' + (i + n A - ^m) f + U - YTY^) f 2 -\( l + x ) Q m = 0, (47) 



where 



X = S = 5 (a = 0)exp fdlna^j 



(48) 



Unlike for ACDM, the higher order terms in S of Eq. (|20| in the Szekeres models force us to keep the integral form 
for X. One could attempt to integrate this numerically as has been done in the past. Alternatively, as in our case, 
we can simply observe that / is related to the growth rate G in Eq. (j2tl)) as 

f = l + a^, (49) 

where G' and G, and thus /, are all functions of a, f2^, and S1 A , and we require = 1 — fi^ — . Thus the function 
/ can be calculated from G and G', which have already been integrated through the Runge-Kutta method used in 
the first sections of this work. In the next subsection, we compare the function / in the Szekeres models directly to 
measured values of /. 

Also important in studies of the growth of structure is the growth index 7. The growth factor / was shown to be 
well approximated by }69h72i | 

/ = (50) 
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FIG. 10: Plots of the best fit growth factor / to observations for Szekeres and ACDM models. Best fit parameter values are 
determined by \ 2 minimization with and without priors. The main result plotted here is that the Szekeres model with no 
priors fits the growth data with \ 2 almost indistinguishable from those of the ACDM. The best fit Szekeres model requires less 
matter density and more spatial curvature. For further exploration, priors on the Szekeres model are applied and consist in 
constraining the total matter density to be that of only baryons, with values as limited by the BBN (i.e. 0.039 < f^" < 0.049). 
Also, the priors when imposed on ACDM consist in fixing the density parameters to those of the concordance model with 
Q.° n = 0.27 and Q.° A = 0.73. The plots are shown with growth factor data both including (left panel) and excluding (right panel) 
available Lyman-a data for the parameter sets described in Table HIl 



More recently, the growth index has been demonstrated to be a useful parametrization of /, capable of distinguishing 
between ACDM and other alternative dark energy and modified gravity models, having different characteristic values 
in each [52, 62-68]. It can be written in ACDM with a dark energy equation of state parametrized by w as an 
expansion to lowest order in 1 — Q m and Qk 64] 

= 3(^-l)(l-n m )-(3^ + l)n fc ( 
7 (6iy-5)(l-fi ro )-2(3w + l)O fc ' [ ' 

For early times (fi TO = 1) and zero curvature, we find 7 = 6/11. This result is nearly independent of the choice of 
and Qfc, with the value of 7 changing only at the percent level. However, it has been shown that the original 
ansatz In / = jhin m must be altered in some cases, such as for a curved space. It is thus not immediately clear that 
we can apply directly a relationship like / = fi^ for some 7 in a general Szekeres space, and the derivation of an 
analytic form for such a 7 is not trivial if possible at all. We will explore whether a 7 can characterize the Szekeres 
models in the following section. 



B. Fitting and results 

We use the measured values of / given in Table Q] to constrain the best fit cosmological parameters f2^, il A , and 
f2? using a \ 2 minimization. We define the x 2 for N data points as 



N 

,2 



i E />■>-&.(■..) , (52) 



where tr(aj) is the absolute error in the measurement at a,-. We construct the \ 2 m this way so as to be able to 
directly compare results that use a different number of data points. The \ 2 1S then one when the average difference 
in / from observations is equal to the average error. 

We first calculate the \ 2 with no priors on the values of the cosmological parameters for both the Szekeres models 
and ACDM, and then impose for the Szekeres models the requirement that 0.039 < < 0.049 from BBN [irj ]. 
assuming that there is no dark matter component (O? = 0). For comparison, we also calculate the % 2 for the 
concordance ACDM cosmology with values fi^ = 0.27 and f29 = 0.73. Finally, we determine the parameters for the 
Szekeres models which best fit the growth in the standard ACDM cosmology. We summarize these results in Table HIl 
and compare the resulting /(a) to / bs(^i) m Fig. 1101 
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FIG. 11: The relative error (/ — fi^)// of the three unique parameter fits in Table HT1 for both Szekeres and ACDM models. 
The relative errors of the three ACDM parameter sets overlap and are nearly zero even at late times, while there is significant 
error in the three Szekeres parameter sets, particularly at late times. 

We find that there is negligible difference in fitting power {x 2 ) between our best fit Szekeres parameters with no 
priors and those for ACDM, both including and excluding Lyman-a data. In general, we find that the Szekeres 
models fit the measured / with comparable f2° to that of the ACDM model, but with a lower matter content and a 
larger spatial curvature. In order to explore the Szekeres fitting result further, we include the priors from BBN that 
0.039 < < 0.049 and set f^ m = 0. We find in this case that the Szekeres models perform much worse in fitting the 
measurements, with almost a factor of three increase in x 2 when excluding the Lyman-a data. This indicates that 
while the Szekeres requires smaller values for the dark matter density, it still needs some in order to have a good fit 
to the growth data. As a simple check, we find that the concordance ACDM parameters are still a good fit to the 
growth data but with non-negligible increases in x 2 when compared to the best fit parameter choices. For further 
comparison to the ACDM model, we find the Szekeres models most closely mimic the growth of the standard ACDM 
model with parameters = (0.11,0.71,0.18) and x 2 = 2-4 x 10~ 5 . 

The information in Table ITT1 is represented visually in Fig. [ljjj where the various / fits are generally very consistent 
and fit well to the data, except for the Szekeres fit with the priors listed above, which clearly fits the data very poorly 
both by visual inspection and x 2 value. We also note that in the right panel, where we exclude Lyman-a data, the 
best fit Szekeres / with no priors is able to fit all data points, while the standard ACDM falls just outside the error 
bar on the data point at a = 0.56 (z — 0.8). 

Our finding of a nonzero curvature using the Szekeres models is consistent with other studies of inhomogeneous 
models and averaging studies, which have shown that using a more general, inhomogeneous cosmology strongly impacts 
the determination of cosmological parameters. For example, by introducing inhomogeneities to explore changes in 
the distance to the cosmic microwave background, [48| found that parameter determinations were strongly affected, 
with a positive curvature and sa 0.8 — 0.9 preferred. This is similar to our determination of the parameters 
for the Szekeres models in the case where we require only baryonic matter to contribute to the matter density and 
find a small positive curvature and large fi^. Similarly, [73| explored the effects of averaging inhomogeneities on 
the measurement of cosmological parameters and found that without prior assumptions on A, a substantial negative 
curvature was favored. This is also consistent with our best fit determination of the Szekeres models, which contains 
a significant negative curvature, as well as lower matter density than required in ACDM. Further, [74] found that the 
averaging will in general affect the Friedmann equation by introducing a curvaturelike term. Thus the measurement 
of an FLRW model with negligible curvature, which is considered to be associated with the averaged background of 
a general Szekeres model, may indicate that a nonzero curvature in the true Szekeres model is appropriate. 

In addition to fitting the growth factor, we also attempt to fit a growth index. To check whether we can apply the 
ansatz In / = 7 In Q m to the Szekeres models, we perform a numerical x 2 minimization to find the best fit 7 for the 
resulting parameters in Table HU As a check that this fitting method performs well, we first apply it to the ACDM 
parameter sets, finding as expected that 7 w 0.55. We find, however, that there does not exist a unique 7 for the 
Szekeres models, but instead that 7 varies with different J7 m , and we find very poor x 2 values for these fits. Further, 
there exists only a small subset of parameter choices for which we can fit a 7 at all, indicating that the usual ansatz 
used for FLRW models fail for the Szekeres models. This is demonstrated in Fig. [Til which shows a plot of the 
relative error (/ — fi^J// of the three unique parameter fits in Table [TTI for both Szekeres and ACDM models. While 
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the relative error in the / = relationship for ACDM overlaps for the three parameter sets and is nearly zero even 
at late times, there is significant error in the three Szekeres parameter sets, particularly at late times showing that 
the usual ansatz cannot be used here and direct comparison to the growth factor data is required. 

We conclude then that the Szekeres models are able to exactly mimic a standard ACDM growth history characterized 
by /, and to fit observed values of the growth factor with comparable accuracy to the ACDM model. The Szekeres 
models require comparable to ACDM, but significantly less dark matter, which is offset by a large negative 
curvature. Our parameter determinations are roughly consistent with previous studies of inhomogeneous models and 
averaging. We also find that the growth index parameterization fails for the Szekeres models, where a single-value 
7 is unable to fit /. However, any conclusions based on comparisons of the Szekeres models to observations, and in 
particular the determination of cosmological parameters, should be made with caution. The measured values of / are 
based on assumptions of ACDM, which do not hold in a more general Szekeres model. The cosmological parameters 
described in the preceding sections are defined in the same way as their ACDM counterparts and should represent 
the same fractional energy contents in the universe. However, the energy densities do not have strictly the same 
connections to physical quantities in both cases. For example, the Hubble parameter in FLRW is proportional to the 
expansion scalar 0, but in Szekeres there is also a dependence on the metric function H (see Sec. W\ . Thus, a more 
careful consideration is required when comparing the Szekeres parameters directly to ACDM. 

VII. CONCLUSION 

The standard approach to studying large-scale cosmological structure evolution is to apply linear perturbations to 
a smooth FLRW metric. However, in order to go beyond the linear perturbation framework, it is worth exploring 
alternative exact frameworks. Such exact frameworks can include linear and nonlinear effects from for example large- 
scale superstructures in the lumpy universe that we observe. In this work, we have studied the growth of large-scale 
structure in the Szekeres inhomogeneous cosmological models in the presence of a cosmological constant and allowing 
for nonzero spatial curvature. Using the Goode and Wainwright formulation of the solution, we have identified a 
density contrast for both classes (5 in class I and 8 in class II) in terms of the metric functions F and H . The F 
function obeys a second-order time differential equation like the density contrast in FLRW, and this allows us to write 
evolution equations for 5 and 5 that have linear terms formally identical to the linearly perturbed FLRW equation. 
However, our equations also contain nonlinear terms. Most importantly, unlike linear perturbations, the equations 
here are derived in an exact framework, meaning S and 8 are not constrained to be smaller than 1. 

To study structure growth in class I, we have constructed (today) a simple nonsymmetric and inhomogeneous 
region that physically represents a central void and neighboring supercluster surrounded by an asymptotically smooth 
background. The model matches to an almost FLRW model with standard ACDM parameters at r = 50 Mpc. By 
numerically solving the growth equation for S, we have evolved the structure backward in time to see the region 
smooth out, i.e. become more homogeneous, with both the overdensity and the void shrinking continuously. The 
class I dependence on the coordinate r of the metric functions allows us to model such a single structure, while for 
class II such an implantation is possible via the metric functions /3± since M is constant and a — a(t). 

The class II models, as we have set them up, are described by their input cosmological parameters today and 
describe exact deviations from an underlying FLRW model (obtained by letting 5 = 0). For class I models, these 
input cosmological parameters today can be understood to apply at some relevant r, and we restrict ourselves to the 
independently evolving spatial section defined by that constant r. For both classes, we find a growth up to several 
times stronger than that of the ACDM model. For example, we find that for flat class II models with and Vt\ 
values comparable to those of ACDM, Szekeres models experience significantly stronger growth than their linearly 
perturbed FLRW counterparts. The flat model with = 0.04 provides a qualitatively consistent growth history 
that, if the matter content is interpreted to be solely due to baryons, appears not to require dark matter in order to 
achieve a similar growth rate to ACDM today. This value of agrees with the value of as predicted by big bang 
nucleosynthesis. Flat and curved class II models with nonzero f2° experience suppression of the growth rate after the 
onset of A dominance, as expected. 

We have also derived expressions for kinematic quantities in both classes, including the shear scalar, expansion 
scalar, and electric part of the Weyl tensor, which represents the tidal gravitational field. We explored the behaviors 
of |<t|/0, p, and \E\/p for the class I model along directions through the supercluster and orthogonal to it, as well 
as for class II models with various combinations of cosmological parameters. These quantities do not experience 
singularities in the class I model, and in class II, the inclusion of A is found to remove divergences that were present 
without it. 

To investigate the ability of Szekeres models to match observations, we derived the evolution equation for the 
growth factor / in class II and also expressed it directly in terms of G and its derivative, the latter being more 
convenient for our purposes. We used a x 2 minimization to obtain best fit parameters from redshift space distortions 
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and Lyman-a forest data for Szekeres and ACDM models both with and without priors. For Szekeres, the parameters 
(fi^, , 0£) = (0.12,0.59,0.29), obtained without priors and excluding Lyman-a, achieve the smallest x 2 - Further, 
the growth factor for these Szekeres parameters is able to fit all the data, while the corresponding ACDM curve misses 
the error bound on one point. With no priors imposed, we find that the Szekeres models achieve the same \ 2 as the 
ACDM model but require different combinations of the cosmological parameters. We have also determined that the 
parameterization of / by the growth index 7 does not apply to the Szekeres models in a simple way. For example, 
the usual fi^i ansatz does not work, and it is necessary to directly fit /. 

With these results, we find that the Szekeres inhomogeneous cosmological models can provide a framework to 
analyze the growth of large-scale structure using exact equations not limited to the linear regime or first-order terms 
in the density contrast. Szekeres models subjected to current growth data can provide competitive fits compared 
to the concordance ACDM perturbative scheme, but they require different values for the matter density and spatial 
curvature. Our results in this work will be useful in developing a complete framework based on inhomogeneous 
cosmological models in which current and future observations can be analyzed and interpreted. 
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Appendix A: Evolution of Szekeres Models in the GW Representation 

We describe the time and space dependencies of Szekeres models in the Goode and Wainwright (GW) representation 
HI]. The GW metric has the form 

ds 2 = -dt 2 + a 2 [e 2 »{dx 2 + dy 2 ) + H 2 W 2 dr 2 ] . (Al) 



a. Time dependence 



The solutions to the equations governing the time evolution of Szekeres models have the same form for both classes. 
When A is nonzero, t he g eneral solution of the Fricdmann equation [Eq. ([5])] involves elliptic functions and can be 
found in, for example, [381 ] - We give here the original forms with A = from the Goode and Wainwright paper, which 
are written in parametric from with parameter 77: 

a = M^, with t-T = Mh, (A2) 

where 

77 — sin 77, k = +1 
h{rj) = { sinh?7-7i, k = -1 (A3) 
ri 3 /6, k = 0. 



r\ varies over < r\ < 2n for k = +1 and < r/ < 00 for k — 0, —1. By Eqs. (| A2|) and (|A3[) . we see that the scaling 
function a has the same time dependence as that of an FLRW dust model. As was pointed out in Sec. UH a condition 
to have a physically meaningful solution is that M > [23J, and when k = 0, — 1 we take a > 0, which corresponds to 
an expanding universe. 

The two independent solutions to the second-order differential equation in F [Eq. (|3J)], are 

(6M/0) [1 - (77/2) cot(r?/2)] - 1, k = +l 
f+ = { (6M/ a )[l-(77/2)coth(77/2)] + l, k = -l (A4) 
r/ 2 /10, k = 
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and 

(6M/0) cot(7y/2), fc = +1 

/_ = { (6M/a) coth(?7/2), fc = -1 (A5) 

24/?7 3 , fc = 0. 

fc. Space dependence 

Unlike the time dependence, the space dependence is different for each class. 

class I: a = a{t,r), a, r ^0, f± = f±{t,r), T = T(r), M = M(r), 

with 

e v = f{r)[a(r)(x 2 + j/ 2 ) + 2b(r)x + 2c(r)y + d(r)] -1 , (A6) 

e/4 = a(r)d(r)-6 2 (r)-c 2 (r), e = 0,±l, 

A = fv, r -k0 +) W 2 = (e-kf)-\ 

8+ = -fc/M, r /(3M), p_=fT, r /(6M), (A7) 

and where a comma denotes differentiation with respect to the variable that follows. The functions f(r), a(r), b(r), 
c(r), d(r), and T(r) are smooth but otherwise arbitrary. The function a(r) is distinct from the generalized scale factor 
a(t,r). 

class II: a — a(t), f± = f±(t), T = const, M = const, 

with 

e u = [l + k/A(x 2 + y 2 )]~ 1 , fc = 0,±l, W = 1, (A8) 

, f e "{ a (r)[l-|(i 2 + y 2 )]+fe(r)£ + C (r)y}-fc/3 + , fc = ±1 . . 

\ a{r)+b(r)x + c{r)y- 0+(x 2 +y 2 )/2, fc = 0. 1 J 

The functions a(r), b(r), c(r), and /3±(r) are smooth but otherwise arbitrary. The function a(r) is distinct from the 
scale factor a(t). 

Appendix B: Relationships Between the GW and LT-like Class I Metric Representations 

The Szekeres metric can be written alternatively in the LT-like representation [75j j 

ds 2 = -At 2 + ($ 7_y )a dr' + fj(dp 2 + V), (Bl) 
where $ = $(t,r) represents an areal radius and is defined by 



(^t) 2 = -k(r)+ 2 -^l + ^ 2 1 (B2) 



where M(r) corresponds to the active gravitational mass within a sphere of coordinate radius r, and 



£{r lP ,q) = i 



S(r) + S(r) +eS[r \ 



(B3) 



The constant e = 0, ±1 describes the geometry of the (p, q) 2-surfaces. 

If we consider a curvature of the form k(r) — Kf 2 (r) (K = \k(r)\/k(r) — 0, ±1) as described by [43[, then we can 
identify f(r) = \/\k(r)\ as the GW metric function f(r). This holds except when k(r) = (K = 0), in which case 
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f(r) 7^ is arbitrary. We can then build a complete correspondence between the metric functions in the following 
two representations: 



a(r) = 
b(r) = 
c(r) = 
d(r) = 



2S(r) 

~P{r) 
2S(r) 

-Q(r) 
2S(r) 

P 2 (r) + Q 2 (r)+eS 2 (r) 



2S(r) 



f(r) 
k 
e 

a(t,r) 

M(r) 
W(r) 

H(r) 



m 

K 
e 

m 

M(r) 

iHr) 
1 



(B4) 



$, r (t,r) 



£,r {r,P,q) 



The above relations have some consequences for the geometrical interpretation of the GW coordinates, following 
insights provided by the LT-like representation. We can immediately identify that the (i , y) and (p, q) 2-surfaces obey 
the same geometry in both metric representations, as indicated by e, and that the GW spatial curvature (k) in some 
region matches that of the LT-like spatial curvature (K = k(r) / \k(r)\) . If we write the GW e v in a suggestive way as 
e " = f(r)/E(r,x,y), where [H 



E(r,x,y) = - 



(x + b/a) 2 (y + c/a) 2 



1 /2a 



l/2a 



(B5) 



we can identify the GW x and y as projected coordinates following the same geometry as the LT-like p and q (see 
IT! . §19.5.3), where the GW functions (b/a, c/a, l/2a) take the place of the LT-like functions (P, Q, S). We can then 



translate the GW metric coordinates (t, r, x, y) to both spherical (i, r, 8, 
Within each surface of constant t and r, these transformations are 



and Cartesian-like (t, X, Y, Z) coordinates. 



{ax + b,ay + c) 
(ax + b,ay + c) 
(ax + b,ay + c) 



and 



2~ cot \y ( cos ^' sin( 
(cos<j), sin</)) 
9 ' 

-cothf - | (('o>, 



X = r sin 6 cos < 
Y = r sin 9 sin <j 
Z = r cos 9. 



(e 



H) 
0) 

-1) 



(B6) 



(B7) 



One can see that with a simple substitution following Eq. (IB4[) . Eq. (|B6[) reduces directly to the transformations for 
the LT-like metric coordinates, given, for example, in |43l |. 

Finally, it is clear that a(t,r) and M(r) play the same physical roles (albeit appropriately rescaled) as their well- 
studied counterparts $(t, r) and M(r) in the LT-like representation. 
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FIG. 12: The function M(r) = f s (r)M(r) compared to the background FLRW M b (r) = p b r A /6. M(r) grows less quickly than 
Mt(r) at r < 15 Mpc due to the central void, but is then compensated by an overdensity and matches to Mj,(r) at r > 50 Mpc. 



Appendix C: Modeling Exact Inhomogeneous Structures in the Class I GW Representation 

For studying the growth of structure in the class I GW representation of the Szekeres metric, we construct the 
metric functions by considering a simple inhomogeneous and nonsymmetric structure similar to that of Model 1 in 
[29l ] . It physically represents a central void (underdensity) with a neighboring supercluster (overdensity) . Constructing 
specific models to represent cosmic structures seems to be more intuitive in the LT-like representation due to a simpler 
comparison to LTB models and some previous work using this LT-like representation has been done [30L HH, [39| . 
Therefore, we perform initial calculations for the model in this representation and then convert the metric functions 
to the GW representation using the relationships described in Appendix B. 

Our model is matched to an almost FLRW background model described by = 0.27, f2° = 0.73, and Ho = 72 
km s _1 Mpc -1 at a radius of 50 Mpc. This background FLRW model has a homogeneous background density given 
by Pb — QmPcr- The model is fully defined through the following series of simple steps: 



(i) We begin by constructing the function 

ni _ / / 9* \ 3 \ 

(CI) 



MM = f(r)M(r) = —r 3 [ 1 - exp 
6 



-3(1 
a 



as shown in Fig. [TJJ where a = 30 Mpc, such that at a radius of 50 Mpc we match to the almost FLRW 
background density. For comparison, M\,{r) — p^r 3 /6 of the FLRW background density is also shown. 

(ii) The metric functions P, Q, and S are chosen. In our case, we choose the functions to match those used in Model 
lof[H|: 

P(r) = 10, 

Q(r) = -1131n(l + r), (C2) 
S(r) = 140. 

These define the function £(r,p 7 q) given in Eq. (|B3|) . which governs the (p,q) dependence of the resulting 
density distribution of the structure, and e = +1 for this model. 

(iii) The r coordinate is defined such that &(t a , r) — r, where t is the age of the universe today. 

(iv) We integrate Eq. (|B2[) in order to find k(r), 

d$ 



I 



^_ fe(r) + ™p + ^ 



= t , (C3) 



where we have chosen the bang time to be everywhere zero (i.e., £s(r) = 0). 
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FIG. 13: The density distribution (p) given by Eq. © of the model described in Appendix C is shown in units of the background 
density (pb)- The density is plotted in Cartesian-like coordinates for the Z — plane [77l ]. converted from the metric coordinates 
through Eqs. ()B6|I and ()B7[) . The model consists of a central void (underdensity) with a neighboring supercluster (overdensity). 



(v) The function $(t,r) is now fully determined by Eq. (|B2|) . 

(vi) Now that M(r), k(r), and $(i, r) are known for the model, we can freely transform to the GW metric represen- 
tation following the relations given in our Appendix B. 

The density distribution of the model is given by Eq. J6j, where S(t,r,x,y) is evaluated as discussed in Sec. IHI Al 
The resulting density distribution in the Z — plane [771 ] is shown in terms of the background density in Fig. 1131 
Though the model is designed to match a ACDM background with zero curvature, the Szekeres region itself has a 
nonzero curvature. 
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